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Abstract. We consider super-Brownian motion whose historical paths reflect from 
each other, unlike those of the usual historical super-Brownian motion. We prove tightness 
for the family of distributions corresponding to a sequence of discrete approximations but 
we leave the problem of uniqueness of the limit open. We prove a few results about path 
behavior for processes under any limit distribution. In particular, we show that for any 
7 > 0, a "typical" increment of a reflecting historical path over a small time interval At is 
not greater than (At) 3 / 4-7 . 

1. Introduction. 

The present article has been inspired by two probabilistic models — superprocesses 
with interactions and reflected particle systems. 

The first person to study a reflecting system of particles was Harris [H] who considered 
an infinite system of Brownian particles on the line. He proved that if the initial positions 
of the particles are points of a Poisson point process, then for a large time t the distribution 
of a single particle is normal with the standard deviation (2t/-7r) 1//4 . Spitzer [S] analyzed 
a similar model with particles moving along straight lines between collisions. See [DGL1, 
DGL2, G, Ho] for related results. 

The simplest superprocesses, for example, super-Brownian motion, are continuum lim- 
its of branching systems in which the branching mechanism is independent of the positions 
of particles. There has been considerable activity studying models with interactions. Many 
articles are devoted to models with catalysts, see, e.g., [DF, De]. Various other models 
with interactions are discussed in [AT, BHM, EP, P3]. See in particular [P4] and references 
therein. 

We will study a model similar to that introduced by Harris, in that we will start with 
linear Brownian motion as the spatial process. We will attempt to build a corresponding 
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superprocess with historical paths that do not cross over although they may touch each 
other. 

Our construction is based on a sequence of discrete approximations. Consider for 
every e G (0, 1] a branching particle system which starts initially with N e particles located 
respectively at x\ < ■ • ■ < x e N ^ . Particles move independently in space according to linear 
Brownian motion and are subject to critical binary branching at rate e~ x . To be specific, 
the lifetimes of the particles are exponential with parameter e~ x and when a particle dies 
it gives rise to or 2 new particles with probability 1/2. 

Let us now introduce our basic assumptions. Let 



H e := e^<5 



and assume that there is a finite measure |jonl such that 

\i e ->/i, (1.1) 

£— >0 

where the notation (w) indicates weak convergence in the space M/(R) of finite measures 
on R. In addition, if supp [i denotes the topological support of [i, we assume that 

supple > supp/U, (1-2) 

in the sense of the Hausdorff metric on compact subsets of R (in particular, we assume 
that supp fx is compact) . 

Let Xf denote the random measure equal to e times the sum of the Dirac point masses 
at the positions of particles alive at time t. Then, 

(X!,t>0)^(X t ,t>0), (1.3) 
e— >o 

where the limit process is super-Brownian motion in R with branching rate 7 = 1 (through- 
out this work we consider only this branching rate) and initial value //, and the convergence 
holds in distribution in the Skorohod space D(R + , M/(R)). The convergence (1.3) is the 
standard approximation of super-Brownian motion (see e.g. [P4]). Note that assumption 
(1.2) is not needed for (1.3) but it guarantees that the graph of X 6 also converges in dis- 
tribution to the graph of X (see Lemma 2.3 below), a property that plays an important 
role in our arguments. 

For each particle alive at time t, we can consider its historical path, which is the 
element of C([0,£],R) obtained by concatenating the trajectories of the ancestors of the 
given particle up to time t. Denote by Yf the historical measure equal to e times the sum 
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of the Dirac point masses at the historical paths of the particles alive at time t (Y t s is thus 
a random measure on the set C([0, t], R) of continuous mappings from [0, t] into R). Then 
the convergence (1.3) can be reinforced as 

(l?,f >0)-^->(Y„f >0), (1.4) 

where the limit process is now historical super-Brownian motion started at \i. 

For every e > 0, we can use the original branching particle system to construct a 
new system with reflection. The branching mechanism (critical binary branching at rate 
e -1 ) is the same as in the original system, but the particle paths in the new system 
reflect against each other. A precise construction is given in Section 3, but let us give 
an informal description. The reflected system is such that for every t > 0, the set of 
positions of particles at time t is the same as in the original system, and in particular the 
branching times are the same. During the time interval between and the first branching 
time, the vector of positions of the particles labeled 1,2, ...,AT e in the reflected system 
is the increasing rearrangement of the vector of positions of the particles in the original 
system. Suppose that at the first branching time, denoted by £, a particle dies and gives 
rise to 2 children. If the location of this particle is the j-ih coordinate in the increasing 
rearrangement of the vector of positions at time £— , we will say that in the reflected 
system particle j has given rise to two children labeled jl and j2. Then on the interval 
between £ and the second branching time, the vector of positions of the particles labeled 
1, . . . , j — 1, jl, j2, j + 1, . . . , N £ in the reflected system is again the increasing rearrangement 
of the vector of positions of the particles in the original system. We can easily continue 
this construction by induction. 

Denote by Xf and Y t e the analogues of Xf and Y t £ for the the system with reflection. 
We have Xf = Xf since the set of positions of particles is the same at every time t in both 
systems. On the other hand, Yf is typically very different from Y t s . Indeed, the following 
property holds for any two paths w, w' in the support of Y t e : Either w(r) < w'(r) for every 
< r < t, or w(r) > w'(r) for every < r < t. 

The main purpose of this work is to try to understand the limiting behavior of the 
branching particle system with reflection as e — > 0. Our primary objective was to get an 
analogue of the convergence (1.4) when the processes Y s are replaced by Y s , giving infor- 
mation about the individual paths in the system with reflection. We did not completely 
succeed in this task, but we can prove the following result, where W denotes the set of all 
stopped paths, or equivalently the union over alH > of the sets C([0, £], R). 
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Theorem 1.1. Let S be a sequence of positive numbers converging to 0. The laws of the 
processes (Y t s ,t > 0) for e G £ are tight in the space of all probability measures on the 
Skorohod space D([0, oo), M/(VV)). Furthermore any limiting distribution is supported on 



Hence, by extracting a subsequence if necessary, we can assume that the sequence of 
processes Y e converges in distribution towards a process Y with continuous paths with 
values in M/(W). Note that, for every t > 0, the measure Y t is supported on C([0, £],M). 
Although the question of uniqueness of the limit remains unsolved, we are able to derive 
several results on the path behavior of the process Y. 

First note that, since Xf = Xf for every t > 0, the convergence (1.1) implies that the 
M/(M)-valued process X defined by 



is a super-Brownian motion started at fi. In particular, it is known (see [KS], [R]) that a.s. 
for every t > the measure X t (dy) has a density denoted by x t (y), and that there exists 
a jointly continuous modification of (xt(y), t > 0, y G M). 

The next result shows that for any 7 > 0, a typical oscillation of a reflecting histor- 

3 

ical path is not greater than (At)* -1 , and hence much smaller than a typical Brownian 
oscillation (At) 2. This result is consistent with the Harris [H] estimate, if we translate the 
large-time asymptotics to small-time asymptotics. 

Theorem 1.2. Almost surely for every t > and every r G (0,t), for every path w G 
supply, the condition x r (w(r)) > implies that, for every 7 > 0, 



A more precise version of Theorem 1.2 is given in Section 5 (Theorem 5.10). It is not 
hard to check that if we fix t > and r G (0, t) (fixing r is in fact enough), the condition 
x r (w(r)) > 0, and thus the conclusion of the theorem, will hold for every path w G supp Y t , 
a.s. Alternatively, for every fixed t > 0, the conclusion of Theorem 1.2 holds for a set of 
values of r G (0,£) of full Lebesgue measure, for every w G supply. We believe that 5% 
is the "typical" size for the oscillation of a historical reflected path although we have no 
lower bound justifying this claim. 

We also study the behavior of reflected historical paths at a branching point. If w and 
w' are two reflected historical paths that coincide up to time r > (meaning informally 
that the corresponding "particles" have the same ancestor up to time r), we show that the 



C([0,oo),M/(>V)). 




lim sup 

540 



w(r + 5) — w(r)\ 



= 0. 
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distance between w(r + 5) and w'(r + S) grows linearly as a function of 5, up to logarithmic 
corrections. The precise statement is as follows. 

Theorem 1.3. Let t > 0. If w and w' are two distinct elements of C([0, t], R), we set 

7 WjW / = inf{r > : iu(r) 7^ w'(r)}. 
Then a.s. for any two distinct paths w,w' G supply such that 7 W)U ,' > 0, we have 



\w{lw,w> + 8) - w'{^ w>wl + J) I _ 
lio* 2<51og|log<J| 



hmsup ^ ^ , I ^ = g (w(7w;W/)) > 



and, for every 7 > 0, 

Tto <5 1 log 5 1 _1 



,. Mlw,w' + $) - W'(lw,w' + $)\ 
lim ; — ; ri — ; ' = OO. 



Our proofs rely on several known results on super-Brownian motion. In particular, 
we use the Brownian snake idea [L2] in an essential way, both in the proofs and for giving 
more precise versions of the results. For instance, as a key step towards Theorem 1.1, we 
get a uniform continuity result (Theorem 4.1) for the historical paths of the approximating 
branching particle systems with reflection. The proof of this result requires some precise 
information about the genealogical structure of the approximating systems, which seems 
to be more easily accessible via the snake approach (cf Lemma 2.1 below). 

For an introduction to the theory of superprocesses (measure- valued diffusions) and 
historical processes, the reader may consult [Da, Dy, DP, L2, P4]. 

The paper is organized as follows. Section 2 describes the specific coding that we use 
to represent the genealogical structure of the approximating branching particle systems. 
This section also contains a few important preliminary results. Section 3 presents the con- 
struction of the systems with reflection. Tightness results are given in Section 4, including 
a more precise form of Theorem 1.1. Section 5 contains the proof of Theorem 1.2, and is 
the most technical part of the paper. Finally, Theorem 1.3 is proved in Section 6. 

We are grateful to Carl Mueller, Ed Perkins, Tokuzo Shiga and Roger Tribe for very 
useful advice. 

2. Coding discrete trees. 

We will describe a method that provides a coding of the genealogy of the branching 
particle systems introduced in Section 1, in a consistent way for all values of the parameter 
e G (0, 1]. This method involves embedding branching trees in a path of reflected Brownian 
motion, and is based on [LI] (see also [NP]). 
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2.1 Markov chains embedded in reflected Brownian motion. 

Let P = (ft, s > 0) be distributed as twice a reflected Brownian motion on R + : 

(ft,s>0) ( = } (2|B a |,s>0), 

where B is a standard linear Brownian motion, with So = 0. The reason for the factor 2 
will be clear later. We denote by (Lg, x > 0, s > 0) the jointly continuous family of local 
times of ft normalized in such a way that, for every nonnegative Borel function ip on R_|_, 



/ (p(P 8 )ds= / (p(x)L^dx. 
Jo Jr+ 



Also set r r = inf{s > : L° s > r}, for every r > 0. 

For every e e (0,1], we introduce a sequence of stopping times (T%,k = 0,1,...) 
defined inductively as follows: 

T £ = inf{s > : ft = 2 £ }, 
T| fc+1 = inf > T| fc : sup ft - ft = 2 £ }, 

T| fc+2 = mf{u > T| fc+1 : ft - inf ft = 2 £ }. 

It is simple to check that the variables Tq,T[ — Tq,T| — Tf , . . . are independent and 
identically distributed. To see this, note that if {"it it > 0) is a reflected Brownian motion 
with initial value 70 = b > 0, the process 

"it - inf 7 S 

0<s<t 

is again a reflected Brownian motion, with initial value 0, and also observe that ftps > 2e 
for every k. 

As E(Tq) = e 2 , standard arguments show that for every K > 



sup 



T f s /e 2 ] - S 



0. 



(2.1) 



(First establish this convergence along the sequence e n = n 2 and then use monotonicity 
arguments.) Thus, 

0. (2.2) 



sup 

s<K 



-ft 



e^0 



For k = 0, 1, . . ., set 



'2fc 



2fc + l 



It is easy to verify that k = 0, 1, 2, . . .) is a time-inhomogeneous Markov chain with 
values in R + , whose law can be described as follows (see [LI] Section 3 for details): Sq = 
and <S| fc+1 has the same distribution as Sf/c + U, where U is an exponential variable with 
mean 2e, independent of <Sffc> nas the same distribution as — V^) + where F 

is exponential with mean 2e, independent of 'S'ffc+i • 

From (2.2), we have a.s. for every K > 0, 



sup 



%0. 



We then define a continuous-time process (/3f , s > 0) by setting 



^ e2 = ^forfe = 0,l,2, 



and by interpolating linearly on intervals of the form [ke 2 , (k + l)e 2 ]. It is obvious that we 
also have 



supl/^-ftl^UO. (2.3) 



2.2 The correspondence between excursions and trees 

With each excursion of f3 e away from 0, we can associate a marked tree representing the 
genealogical structure of a Galton- Watson branching process with critical binary branching 
at rate e , starting with one individual (the ancestor) at time 0. Here a marked tree 
consists of the set T of edges (i.e., particles), which is a subset of 



U := [J {1, 2} n (by convention, {1, 2}° = {0}), 

n=0 



and the family (£ u ,u E T) of lengths of edges (i.e., lifetimes of particles). 
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Figure 1. 



This correspondence is explained in Fig. 1 for the first excursion of (3 £ away from 0. 
Informally, if (ie 2 ,je 2 ) is the interval corresponding to an excursion of /3 s , the lifetime 1% of 
the individual at the root of the associated tree is the minimum of f3 e over [(i+l)e 2 , (J— l)e 2 } 
and this individual has two children if and only if j — 1 > i + In that case, by decomposing 
the excursion restricted to [(i + l)e 2 , (j — l)e 2 } at the time of its minimum over this interval, 
we get two new excursions, each of which codes the genealogical structure of descendants of 
one of the ancestor's children. The construction of the tree is then completed by induction. 
Note that each time of the form ke~ 2 in the interval (ie 2 : je 2 ) corresponds to one edge of 
the tree (for instance the time of the minimum over [{i + l)e 2 , (j — l)s 2 } corresponds to 
0, see Fig. 1). We refer to [LI] Section 2 for a more precise description and a proof that 
this construction yields the family tree of a Galton- Watson branching process with critical 
binary branching at rate e -1 . (We can now explain the factor 2 in the definition of (3: We 
want the branching rate to be e~ x and not (e/2) -1 .) 

There is a one-to-one correspondence between excursions of /3 £ away from and ex- 
cursions of f3 away from with height greater than 2e: If ke 2 is the beginning of an 
excursion of /3 £ , then T| is the hitting time of 2e by the corresponding excursion of /3. As 
in Section 1, consider for every e G (0, 1] an integer iV e > 1 and assume that the family 
(eN £ ,e e (0, 1]) is bounded and that sN e converges to a > as e — >• (this follows from 
(1.1) with a = (/i, 1)). Let r e denote the A^ £ -th return of f3 e to 0. From the previous 
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observations, (2.1) and the standard approximation of Brownian local times by upcrossing 
numbers, we have 

lim t s = r a , a.s. 

We will write r = r a for simplicity. 

On the time interval [0, r e ], the process (3 £ makes N e independent excursions away 
from 0. These excursions can be viewed as representing the genealogical structure of the 
branching particle system introduced in Section 1. The set of edges, denoted by T e , is then 
a random subset of {1, . . . , N e } x U and conditionally on T e , the corresponding lengths are 
independent exponentials with mean e. The function , s G [0, r e ]) can be reconstructed 
from this collection of marked trees as shown by Fig. 1. Notice that for this reconstruction 
to be possible, it is essential to order the trees and the different edges of every single tree. 

2.3 Discrete and continuous local times 

One reason for considering the processes f3 £ comes from their relation with the up- 
crossing numbers of (3. We first define the (discrete) local times of [3 s . For every x > 
and s > 0, we define 

L e s ,x = e Cardjr G [0, s) : (3 e r = x and (3 e u > x for u G (r, r + 5], for some 5 > 0}. 

In other words, £ _1 Lf' x is the number of upcrossings of f3 e above level x before time s. 

Let Mg(x) denote the number of upcrossings of (3 from x to x + 2e completed before 
time s. More precisely, M|(x) is the number of pairs (w, v) with < u < v < s, such that 
(3 U = x, (3 r > x for every r G (w, v) and v = inf{r > u : j3 r > x + 2e}. 

Then, a.s. for every x > and every integer k > 1, we have 

This identity is easily verified by induction on k (the sequence of stopping times (T|) was 
designed for this property to hold). See also Proposition 7 of [LI]. 

Lemma 2.1. We have a.s. 

lim ( sup sup \L e 8 % T e - L* AT \) = 0. 

v s>0 x>0 J 



Proof. We first observe that a.s. 

lim (sup sup \eM £ sAT£ (x) - L x sAT \) = 0. (2.5) 

v s>0 x>0 7 
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For a fixed value of x, this is nothing but the classical approximation of Brownian local 
time by upcrossing numbers, and excursion theory provides precise estimates for the rate 
of convergence. Using these estimates and monotonicity properties, it is then an easy task 
to prove (2.5), i.e., the uniform version of the claim. 

The statement of the lemma is now a simple consequence of (2.1), (2.4) and (2.5). □ 

Remark. As an immediate consequence of Lemma 2.1 and the joint continuity of Brownian 
local times, we get that 

lim ( sup sup \L e ;* T e - ^aUI) = 0, a.s. 

\x-x'\<5 

Later, we will consider for every e G (0, 1] a process f3 £ with the same distribution as f3 £ . 
If L £,x denote the discrete local times of /3 e , the last convergence still holds in probability 
when L £ s ' x is replaced by L £ 8 ' x (and t £ by r e , with an obvious notation). 

2.4 Branching particle systems and discrete snakes. 

We now consider the branching particle system of Section 1, starting with iV e particles 
located respectively at x\ , x 6 ,, . . . , x £ N ^. We may and will assume that the genealogy of the 
descendants of particle k (present at x £ k at time 0) is given by the tree associated with the 
/c-th excursion of f3 £ (cf subsection 2.2). We will refer to this system as the e-system of 
branching Brownian motions. 

For our purposes, it will be convenient to view the collection of paths traced by the 
branching particles as the range of a path-valued process called the discrete snake. 

By definition, a stopped path is a continuous mapping w : [0, Q — > R, where ( = 
( w > is called the "lifetime" of w (it is convenient to talk about the "lifetime" of a path 
although for technical reasons the path is stopped rather than killed). Let W be the set 
of all stopped paths. Then W is a separable complete metric space for the distance 

d(w, w') = \( w - ( w >\ + sup \w(t A ( w ) - w'(t A ( w >)\. 

t>o 

For any x G R, we write x for the trivial path such that ( x — and x(0) = x. 

With every s G [0, r e ] we now associate a stopped path W £ G W with lifetime (3 £ s . If 
s G [0, r e ) fl e 2 N and (3 £ s = 0, then s is the starting time of the /c-th excursion of f3 £ away 
from 0, for some k G {1, . . . , N e }. We then set W £ = x%. For definiteness, we also set 
W £ e = x £ N ^. If s G [0, t £ ) fle 2 N but (3 s > 0, we can associate with s a unique edge of the k- 
th tree, k being the number of the excursion straddling s. We then let W £ be the historical 
path of the particle in the system of branching Brownian motions that corresponds to this 
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edge. Notice that the death time of this particle is (3 s , and thus (w? = fit- Finally if 
s G [0, t £ ] but s £ £ 2 N, we find an integer j such that je 2 < s < (j + l)e 2 , and let I = j if 
P £ J£ 2 > P £ (j+1)£ 2, but / = j + liip £ j£2 < P £ {j+1)e 2. Then we let W £ be the path Wfa stopped 
at time (3 £ s . 

It is easy to see that conditionally on ((3 £ a , s > 0) the process (W ke2 , < /c < t £ /e 2 ) 
is Markovian. To describe its conditional distribution, let A; G {0, ... , t £ /e 2 } and suppose 
that P £ k+r)e 2 > (otherwise W £ k+Vje2 = x|, if (fc + l)e 2 is the starting point of the j-th 
excursion of f3 e ). If /?( fc+1 ) e 2 < fi e ke 2 (which occurs if k is odd) then W £ k+1 ^ £2 is simply the 
restriction of W^ g2 to [0, /3^ fc+1)e2 ]. On the other hand, if /^ fc+1)e 2 > P £ ks 2, then W r ( e fc+1)e2 is 
obtained from W kg2 by "adding at the tip of W kg2 " a Brownian path of length P £ k+1 ^ £ 2 —file 2 
independent of (W £ £2 ,j <k). 

The following snake property is a consequence of the definition of W £ : If s < s' and s 
and s' belong to the same (open) excursion interval of f3 £ away from 0, then W £ (t) = W £ ,(t) 
for every t G [0, mf ue[SjS ,] (3 £ u \. 

2.5 Convergence to super-Brownian motion 

As in Section 1, we let be e times the sum of the point masses at the positions of 
the particles alive at time t in the e-system. This is equivalent to writing 

X t = / dL 6 / S W e( t y 

Jo 

To justify this formula, recall the correspondence between excursions and trees described 
in Subsection 2.2 and note that each upcrossing time s of /3 e above level t corresponds to 
one particle alive at time t, whose position is W £ (i). Similarly, the historical process Y t £ is 

Y £ = [ T£ dL £ /S W! . 
Jo 

Recall our assumptions (1.1) and (1.2) and the convergence result in (1.3). We next 
prove a result about the uniform modulus of continuity for the paths W £ . For convenience, 
we make the convention that W £ (t) = W £ {j3 £ s ) when t > (3 £ . 

Lemma 2.2. Let r\ G (0, \ ). Then, 
hm( iirf P[\W £ (t + r) - W £ (t)\ < r^~ v , for every t > 0, r G [0,5], s G [0,r e ]]) = 1. 

Remark. This is of course reminiscent of the uniform modulus of continuity for historical 
paths of super-Brownian motion. This lemma is therefore very close to the results of [DIP] 
and [DP], which however use different approximations. 
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Proof. Obviously it is enough to treat the case when x\ = • • • = x e N = for every 
s. We then use an embedding technique that will also play an important role later. Let 
(W S7 s > 0) be the Brownian snake of [L2] driven by the process (/? s ,s > 0) and with 
starting point 0. Recall that this is a continuous Markov process with values in Wo := 
{w G W : w(0) = 0}, whose law is characterized by the following properties: 

• For every s > 0, the path W s has lifetime (3 8 . 

• Conditionally on (/? s , s > 0), the process (W s , s > 0) is time-inhomogeneous Markov, 
and its transition kernels are characterized as follows. If s < s', we have W s >(t) = W s (t) 
for every t < m(s,s') := inf[ S)S /j f3 r , and (W s '(m(s,s / ) + r) — W s > (m(s, s')), < r < 
j3 s i — m(s, s')) is a Brownian path independent of W s . 

Now, for every e £ (0, 1], we may assume that the spatial motions of the particles are 
chosen in such a way that, for every e > and every fc G {0,1,..., r e /e 2 }, 

WL 2 = W T c if k is odd, 

(2.6) 

W^ S 2 = W T * I [0, Pre - 2e] if is even, 

where the notation Wt* \ [0, (3t* — 2s] means that the path ^ is restricted to the interval 
[0, Pt^ — 2e] = [0, /5| e2 ]- In fact, it is immediate to verify that the process (W£ e2 , < /c < 
t £ /e 2 ) defined by (2.6) has (conditionally on f3 £ ) the distribution described at the end of 
Subsection 2.4. 

Note that the family (r e ,£ G (0,1]) is bounded a.s. Then the proof of Lemma 2.2 
reduces to checking that, for every K > 0, 

\imP[\W s (t + r) - W s (t)\ < r*"", for every t > 0,r G [0,<5],s G [0,if]l = 1. (2.7) 

S.[0 

This can be easily done using Borel-Cantelli type arguments. Alternatively, we may also 
use the relations between super-Brownian motion and the Brownian snake [L2], and the 
uniform modulus of continuity of [DP]. □ 

The graph Q £ of the £-system of branching particles is defined by 

G £ = cl( (JO} x suppXf) = {W £ {t) : s G [0,r e ],0 < t < (3 £ s }. 
t>o 

We are interested in weak convergence of Q £ towards the graph Q of X, which we define as 

Q = cl (UW xsuppAV). 

We view both Q £ and Q as random elements of the space of all compact subsets of M + x M, 
which is equipped with the Hausdorff metric. 
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Lemma 2.3. We have the joint convergence 

{(xf,t>o),g s ) ((x u t>o),g). 

Proof. We first consider the case when x\ = ■ • ■ = x s N = for every e. Then we can 
suppose that the processes (W/, s G [0, r e ]) are constructed via the embedding technique 
described in the preceding proof. From (2.1) and (2.6), we get 

(W^, s > 0) (WW, s > 0) (2.8) 

e— >0 

in the sense of uniform convergence. Using Lemma 2.1, we get 

x t = / dLl^SwHt) — — — >■ / dL t s 8 Ws ^ t ) = X t 
Jo Jo 

uniformly in t. (The formula for X t is the Brownian snake representation of super-Brownian 
motion, see [L2].) Furthermore, the convergence (2.8) also implies that 

Q £ = {W!(t) :s<r £ ,t<P e s } {W s (t) :s<r,t<[3 s }, 

and the limit is easily identified with the graph Q of X. Therefore we get the statement 
of the lemma in the special case x\ = ■ ■ ■ = x £ N ^ = 0. 

Before proceeding to the general case, let us make one more observation. Fix 5 > and 
write (Vg, s > 0) for a process distributed as an excursion of W e away from conditioned 
to have height greater than 5. (Alternatively, (V/, s > 0) codes the historical paths of the 
e-system starting with one particle at the origin and conditioned to be non-extinct at time 
5.) It follows from the convergence (2.8) that we have also 

(^, s >o)-^(y s ,,>o), 

e— >0 

where the limiting process is an excursion of W conditioned to have height greater than 8. 
As in the first part of the proof, it follows that the graphs of V 6 (defined analogously to 
Q £ ) also converge in distribution towards the graph of V. Furthermore, this convergence 
holds jointly with that of the measure- valued processes associated with V 6 in the same 
way as Xf was associated with W e . 

Let us consider now the general case. Because of Lemma 2.2 and assumption (1.2), 
it is enough to prove that for any fixed 8 > 0, Q e D {[8, oo) x R) converges in distribution 
to Q fl ([5, oo) x M) (and that this convergence holds jointly with that of X 6 ). Let A e 
stand for the set of indices j G {1, . . . , iV e } such that the j-th excursion of f3 £ has a height 
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greater than 5. Note that the events {j E A £ } are independent with the same probability 
2e(2e + It follows that the random measure 

E ^ 

jeA £ 

converges weakly to a Poisson measure with intensity Note that, conditionally on A £: 
Q e fl [5, oo) x R has the same distribution as 

U (((o,^)+yn[^)xi) 

jeA e 

where Qj yE are independent copies of the graph of V s . If follows that the random sets 
Q e fl [5, oo) x R converge in distribution to 

U (((o,^) + ^ -))n[5,oo) xi), 

jeJ 

where ^2j e jS Xj is a Poisson point measure on R with intensity and, conditionally on 
this random measure, the random sets G(j) are independent and distributed according to 
the law of the graph of V. The canonical representation of superpro cesses allows us to 
identify this limiting distribution with that of Q D ([5, oo) x R). Furthermore, using the 
joint convergence of (V s , X s ), it is easy to verify that the convergence holds jointly with 
that of X e . □ 

3. Branching particle systems with reflection 

3.1 Reflection for deterministic paths 

The purpose of this section is to explain, first in a deterministic setting, the con- 
struction of reflected systems. We consider a deterministic branching particle system in 
M analogous to the ones considered above. At time 0, we have N particles located at 
x\, . . . , xn- Each particle moves in R and gives birth at its death to or 2 new particles. 
As in Subsection 2.2, denote by T the genealogical forest of the population, which is a 
subset of {1,...,AT} x U. Each element v = (k,u) in T corresponds to a particle with 
birth time £ v and death time ( v (as in Section 2, we could alternatively consider the life 
durations £ v := ( v — £ v but in this subsection and the next one it is more convenient to 
deal with the birth and death times). The spatial motion of v is a continuous function 
fv '■ Cv] — ► ^ and fv'(^v') = fv(Cv) if v' is a child of v (then £ v > = ( v ). The historical 
path of v is the continuous function w v : [0, Cv] — > R such that, for every t G [0, £„), w v (t) 
is the position at time t of the ancestor of v alive at that time. 
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We assume that the death times ( v , v E T are all distinct, that the system becomes 
extinct after a finite number of generations and that when a particle dies there is no other 
particle at the same location: For every v E T, f v (( v ) ^ fv'(Cv) for every v' E T such that 
€v' < Cv < Cv'- 

We turn to the construction of the reflected system. This system is such that the 
number and positions of the particles alive at every time t are the same as in the original 
system (thus each death time for the reflected system is also a death time for the reflected 
system). However the genealogical forest T will be different, as will be the spatial motions 
f u , u E T or the birth and death times £ u , ( u , uET. 

Set Rq = and denote by R\ < R2 < ■ ■ ■ < Rm the successive death times in the 
original system. For every k E {1, . . . , M}, let T^) be the set of (labels of) particles that 
are alive on the interval [Rk-i, Rk)- We use induction on k to define sets Tik)-, which will 
represent the particles alive on the interval [Rk-i, Rk) in the reflected system, and the 
corresponding spatial motions. 

To begin with, we have = {1, . . . , N}, and we define fj(t) for every t E [0, R±] 
and every j E by requiring (/i(t), . . . , /at(^)) to be the increasing rearrangement of 
(fi(t), . . . , /iv(t)). Note that the mappings /1, . . . , /jv are continuous. 

Suppose that for some k E {1, . . . , M — 1}, we have defined T^) and the corresponding 
paths (f u (t),t E [Rk-i, Rk]), for u E T(k)i m such a way that Card7( fc ) = Card7( fc ), and, 
for every t E [R k _ u R k ]\ 

• The mapping 3 u — > f u (t) is increasing with respect to the lexicographical order 
on f {k) . 

• The values of f u (t) for u E (counted with their multiplicities) are the same as 
those of f u (t) for u E T{k)- 

By definition, one of the particles in T(k)i say u^, dies at time R k - Then there is 
exactly one u (fc ) E T^ k) such that fu (k) (R k ) = f U(k) (R k )- We set 

T(k+i) = {T{k)\{u(k)}) U{u( fc )l,U( fc )2} 
if has two children in the original system, and 

T(k+i) = T (k)\{u(k)} 

if not. Furthermore, let u\ +1 , . . . , u^ 1 1 be the elements of T( k +i) listed in lexicographical 
order. We define f u (t) for every t E [R k ,R k +i] and every u E T( k +i) by requiring that 
(f u k+i (t), . . . , / fe+i (t)) is the increasing rearrangement of (f u (t), u E T(k+i))- Notice that 

when u E fl T{k+i) the definition of f u (R k ) is consistent with the previous step. 
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Finally, the genealogical forest of the reflected system is 

M 

k=i 

The birth and death times £ u , ( u as well as the (continuous) spatial motions f u in the 
reflected system are defined by the requirement of consistency with the construction of 
7( fc )'s. Note the two fundamental properties: 

• At each time t > 0, the positions of the particles (counted with their multiplicities) 
are the same in the original and the reflected system. 

• If u, v G T with u -< v (-< denotes the lexicographical order) then f u (t) < f v (t) for 
every t G [f u ,Cu] n [£,,£«]• 

Historical paths w M , « 6 Tfor the reflected system are defined in a way analogous to 
the original one. If u, v G T and w -< v then ■u; tt (t) < w v (t) for every t G [0, ( u A £„]. 

3.2 A technical lemma 

Let M G {1, • • • , A}, and consider a branching system consisting only of the particles 
labeled 1, . . . , M at time and their descendants. The new genealogy is described by the 
forest 

T' :=Tn({l,...,M} x U). 

From this new branching particle system, we can construct a reflected system by the 
procedure described in Subsection 3.1. We denote by T' the genealogical forest for this 
new reflected system, and by w' v , v G T' the associated historical paths. In general, the 
historical paths w' v will be very different from those obtained by reflecting the original 
system. Under special assumptions however, we can say that some of the paths w' v will 
also be (reflected) historical paths in the original system. 

Lemma 3.1. Let t > and let I be a bounded interval in R. Suppose that w v (r) ^ I 
for every v G T\T' and r G [0, t] . If v G T is such that £ v > t and w v (r) G / for every 
r G [0,t], then there exists v' G T such that (' v , > t andw' v ,(r) = w v (r) for every r G [0,t]. 
The converse also holds: If v' G T' is such that (' v , > t and w' v ,(r) G / for every r G [0, t], 
then there exists v G T such that ( v > t and w v (r) = w' v ,(r) for every r G [0,t]. 

In other words, the first assertion means that the path w v , or rather its restriction to 
[0, t], will still be a historical path for the new reflected system. We leave an easy proof of 
the lemma to the reader. 

3.3 Reflected branching particle systems 
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For every e G (0, 1], we can apply the construction of Subsection 3.1 to the £-system 
of branching Brownian motions. Note that the assumptions that we imposed on the de- 
terministic system hold with probability one for this random system. We write T e for the 
genealogical forest of the £-system, and (£^,u G %) for the lifetimes of particles. The 
notation T e and (£^, u G T £ ) has a similar meaning for the corresponding reflected system, 
which we call the ^-reflected system. Observe that (T £ , (l e u , u G T e )) and (T e ,(^,ueT £ )) 
have the same distribution. This is so because the spatial motions and branching struc- 
ture for the e-system of branching Brownian motions are independent (a tedious rigorous 
justification could be given, but we feel that the result is sufficiently obvious to allow us 
to omit it). Furthermore, Card T e = Card T s . 

We noticed at the end of Subsection 2.2 that the process , s G [0, r e ]) can be 
reconstructed as a measurable function of the marked trees {^T e ,{i%.,u G %)). Hence, 
we can also code the branching structure of the e-reflected system by a random process 
(flg, s G [0, r e ]) which has the same distribution as (f3 £ s , s G [0, r e ]). The fact that Card T e = 
Card % implies that the time r e is also the end of the N s -th excursion of f3 e away from 
0, and thus r e = r e . The discrete local times of (3 s (cf. Subsection 2.3) are denoted by 
(L s s > x ,xe R+, s G [0,r e ]). 

Finally, we can code the historical paths of the ^-reflected system by a discrete snake 
(W £ , s G [0,r e ]) in a way analogous to what we did in Subsection 2.4. Recall that we 
assume x\ < • • • < x e N ^. As in Section 2, if s G £ 2 Nfl [0, r e ) and (3 e a = 0, we set W% = x £ k if 
s is the beginning of the k-th excursion of f3 e away from (and W £ g = x e Ns ). Otherwise, 
if s G e 2 N n [0, r e ) and f3 e s > 0, then (s, f3 e s ) can be associated with a unique edge u of the 
forest T £ , and we let W £ be equal to w e ul the historical path of u. If s G" £ 2 N, we use the 
same interpolation as in Section 2. A fundamentally important property of the process 
(Wg, s G [0, r e ]), from the point of view of our project, is that for s < s', 

W!(t)< W s At), VtG[0,^A^,]. (3.1) 

This follows from our construction and the end of Subsection 3.1. As in the case of Wg, 
we see that if s < s', then 

W!(t) = W! > (t), Vfe[0, mf M 

u€[s,s'\ 

Because at every time the locations of particles are the same in the reflected system 
and in the original one, the random measure 
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coincides with X^. Things are however very different for the historical measure 

Yt = rdLf8~ e . 
Jo 

4. Tightness of the reflected system 

4.1 Uniform continuity of the reflected paths 

Our first goal is to derive an important uniform continuity property for the individual 
paths of the ^-reflected system (Theorem 4.1 below). From the intuitive point of view, 
reflected paths should have smaller oscillations than "free" paths and so this property 
seems to be a straightforward consequence of Lemma 2.2. However the intuition about the 
relationship between moduli of continuity of free and reflected paths is only correct as long 
as we do not have any deaths. To be specific consider p paths wm, . . . ,W(n) all defined 
on the time interval [0, 1], and let W(i), ■ ■ ■ >^U(p) be the corresponding system of reflected 
paths. Then, if we assume that \w^(t) — wu)(t')\ < (p(\t — t'\) for every i = 1, . . . ,p and 
t, t' G [0, 1] and for some nondecreasing function tp, an easy argument shows that the same 
bound holds when the paths wu) are replaced by W(iy 

It turns out that a similar assertion about moduli of continuity is false if paths may 
have different lifetimes. Fig. 2 shows a system of two paths. In the original system, the 
oscillations of paths over the intervals where they are defined are equal to Z\ —y\ and Z2—yi- 
One of the paths in the reflected system goes from y\ to Z2 and so has an oscillation larger 
than the oscillations of the original paths. In this article we consider Brownian particles 
which die at different times so we cannot use known estimates for the modulus of continuity 
of the original (non-reflecting) historical paths in a direct way. We will use them later in 
a different but quite essential way. 




Figure 2. 
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Recall our notation w^, u G T £ , for the historical paths of the ^-reflected system. By 
convention, w e u (t) = w £ u {Q) if t > 

Theorem 4.1. For every r\ > 0, 



lim limsup P 

6^0 V £ ^ 



sup sup \w £ u (t) - w e u (t')\ > rj ) = 



t,t'>0 
\t-t'\<S 



Proof. Let (<^( p ),£( p )) be a sequence in (0, l] 2 converging to 0. We will prove that there 
exists a subsequence (^( p )? £ ( p )) such that: 

lim ( sup sup \wu (p) (t) -Wu (p) (t')\) = (4.1) 

|*-*'|<<5; P) %> 

in probability. Clearly, the statement of Theorem 4.1 is a consequence of this fact. 

We first explain how we choose the sequence ($( p y £ { p ))- By Lemma 2.3 and the 
Skorohod representation theorem ([EK] Theorem 3.1.8), we may, for every p > 1, replace 
the pair (X e M, G £(p) ) by a new pair with the same distribution (for which we keep the 
same notation), in such a way that 

where X is a super-Brownian motion started at ji and Q denotes its graph. Note that 
the genealogical forest T £(p) , the process /3 £ <p) and the historical paths Wu P \ u G % {p) , 
are reconstructed as measurable functions of the new process X^ £ p\ and that it suffices to 
prove (4.1) for the new historical paths. As a consequence of the remark following Lemma 
2.1, we have 

lim ( sup sup \L £ s ^ (p) - 1) = (4.2) 

p^oo V s > t ,t'>0 SAT SAT ' 

\t-t'\<5 (p) 

in probability. We choose the subsequence (^( p ) 7 £ ( p )) so that the convergence (4.2) holds 
almost surely along this subsequence. 

We will argue by contradiction to prove (4.1). If (4.1) does not hold, then on a set A 
of positive probability, we can find a number rj > and a (random) subsequence p k | oo 
such that, if e k := e' (pfc) and S k := ^ pfe) , 

sup sup \w £ u k (t)-w £ u k (t')\>rj. (4.3) 

|t-t'|<«5 fc 
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From now on until the end of the proof, we will assume that the event A holds. By (4.3), 
for every k > 1, there exist u k G T £k , t k , t' k > with \t k — t' k \ < S k , such that 

\ti £ u k k (t k ) ~ ti £ u k k (t' k )\ > V- 

Clearly, we can assume that t k <t' k < 

Recall that the graphs Qe (p) converge to Q in the Hausdorff metric. In particular, the 
set of all pairs (£&, w^ k (tk)) and (t' k , w^ k (t' k )) is relatively compact. By passing to a subse- 
quence, if necessary, we may assume that t k , t' k — > too, (t k ) — > x\ and w^ k (t' k ) — > x^ 
as k — > oo. We have \x\ — x%\ > Tj, and we take x^ > x\ for definiteness. 

We also know that X^ (p) = X^ (p) converges to X t a.s. as p — > oo, uniformly on 
compact subsets of R + . Hence, both sequences X\ k and X e t k converge to X toa: and 

k 

limsupX £k ((-oc,w £ u k k (t k )}) <X too ((-oo, Xl ]), 
hminfX t ?((-oo,^(4)))>X, oo ((-oo,x 2 )). 

We claim that 

liminf (X £ t k ((-oo, < fe fe (t fc )]) - X«* ((-oo, w £ u k k (t' k )))) > 0. (4.5) 

To see this, we use the discrete snake representation of Subsection 3.3. Write s k G [0, r £fe ) fl 
e 2 N for the time associated with the edge u k of T ek in this representation. By construction, 
™ e u k = w ! k : and (3- 1 ) implies 

KH(-oo,wii(t k )]) = £ k dir k \^ Htk) ^ {tk)] > itf k . 

Similarly, we get 

X s /((-oo,wll(t' k )))<Llf\ 

k 

Hence, 

X£((-oo, - ^((-oo, tSSt (**))) > L £ sf k -Ls k / k - (4-6) 

On the other hand, 

^'*' fe |< sup sup |Lf fc >*-L^>*'|. 
se[o,r £k } t,t>>o 

\t-t'\<5 k 

By the convergence in (4.2), which holds a.s. along the subsequence (^( p )> £ ( p ))) the right 
hand side tends to as k — > oo. This and (4.6) give the claim (4.5). 



T £k,tk 
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From (4.5) and (4.4), we get X too ((—oo,xi]) > X too ((— oo, x 2 )) and thus (recall that 
Xi < x 2 ), X tao ((xi,X2)) = 0. This a priori does not imply that {too} x (^1,^2) H £ = 
as there could be a "local extinction" of X at time too in (xi, £2)- However, by Theorem 
1.4 of Perkins [PI], there can be at most one local extinction at a given time, so we can 
choose x\ and x' 2 with x\ < x\ < x' 2 < x 2 such that {too} x x 2 ] H C? = 0. Since £ is 
closed, we have also [too — 5, too + <5] x [^'i, x' 2 ] H = for 5 > sufficiently small. However, 
by construction, for sufficiently large the paths wff and thus also the graph Q £k must 
intersect [too — 5, too + 5] x [^'1,^2]- This gives a contradiction since we know that Q £k 
converge to Q. This contradiction completes the proof of Theorem 4.1. □ 

4.2 Tightness of reflected discrete snakes 

From now on, we restrict our attention to values of e belonging to a fixed sequence S 
decreasing to 0. For convenience, we extend the definition of the discrete snakes W s by 
taking Wg = W^ e = x £ Ne (and thus (3 e s = 0) for s > r e . 

Proposition 4.2. The laws of the processes W e , e G £, are tight in the space of all 
probability measures on D([0, 00), W). Furthermore, if(W s ,s > 0) is a weak limit point of 
this sequence of processes, we have the following properties. 

(i) If p s := , the process (fi s , s > 0) has the same distribution as ((3 sAt , s > 0). 

(ii) Almost surely for every s < s' we have W s (t) < W s r(t) for every t G [0,/3 s A j3 s >]. 

(iii) The set of discontinuities of the mapping s — > W s is contained in the zero set of 
13. Furthermore, if s < s' belong to the same connected component of the complement of 
the zero set, we have 

W a {t) = W s >(t) for every t G [0, inf j3 r ] . 

r€[s,s'] 

Proof. The hard part of the proof is to show tightness. To this end we rely on the classical 
criteria (see e.g. Corollary 3.7.4 of [EK]). We first observe that the compact containment 
condition is a straightforward consequence of Theorem 4.1. In fact, if n > is fixed, then 
for every integer p > 1, Theorem 4.1 and the construction of the discrete snake W e allow 
us to find 5 P > such that, for e G £ small enough, 

sup sup \W e s (t) - W e s (t') \ > 2~ p ] <i 1 2~ p - 1 . (4.7) 

s>0 t,t'>0 -I 
\t-t'<S p 

(Here and later, we make the convention that Wg(t) = VF s £ (/3f) for t > (3 e s .) It is easy to 
see that an even stronger assertion holds, namely, (4.7) is true for all e G £; this can be 
achieved by taking 5 P even smaller if necessary — note that for any fixed value of e we need 
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only consider a finite number of historical paths. Then let H be a compact subset of R + 
containing supple for e G £, and let A > be a constant. The set 

K :={w G W : w(0) eH,(w<A, 

and \w(t) — w(t')\ < 2~ p for every t, t' G [0, ( w ] with \t — t'\ < 5 P and every p > 1} 

is compact, and it follows from (4.7) that 

P[W* £ K for some s > 0] < rj 

provided that A is chosen large enough. 

Recall the definition of the distance d from Subsection 2.4. We set 



6(e,S) = m£ {sup sup d(W^ W^)} , 



( s i) i S,s'e[Si_i,Si) 

where the infimum is over all finite sequences = so < s\ < ■ ■ ■ < s m _i < r e < s m such 
that inf{|sj — Sj_i|; 1 < z < m} > 5}. As a direct application of Corollary 3.7.4 in [EK], 
the proof of tightness will be complete if we can verify that, for every rj > 0, we can choose 
S > sufficiently small so that 

limsupP[0(e, 8) > rj] <rj. (4.8) 

EBe^O 

We now fix rj > and proceed to the proof of (4.8). As a consequence of Theorem 
4.1, we can choose p G (0, rj/5) so small that, for every e G £, 

sup sup < ?1 >!-?■ ( 4 - 9 ) 

s>0 t,t'>0 J 

I*— *'l<p 

Then, by the tightness of the laws of f3 £ (cf (2.3)), we can choose k > small enough so 
that, for every e G £, 

sup |g-g # |<pl >i-5. (4.10) 

s,s>>0 J 

|s — S 1 I <K 

We denote by E £ the intersection of the events considered in (4.9) and (4.10), so that the 
probability of the complement of E e is bounded above by 2rj/5. 

Set 7 = 77/5. Since \i is a finite measure with compact support, we can easily find an 
integer M 1 and a finite sequence of reals y\ < z\ < yi < Z2 < • • • < ?/m 7 < ^m 7) such that: 

• Zi - yi < 7 for every i = 1, . . . , M 7 , 



[j/j, Zj) contains a neighborhood of supp 



i=l 
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• Kivi}) = Ki z i}) = o» and K[y» z i)) > for ever y * = i,---M 7 . 

By the last condition, a 7 := inf{/z([yi, ^)) , i = 1, . . .M 7 } > 0. Furthermore, if £ is small 
enough, 

SUpp^g C [Jbi,^) 
i=l 

and 

Card{j : ^ G [yi,*)} > |j > 1, 

for every i = 1, . . . , M 7 . From now on, we assume that £ G £ is small enough so that the 
last two conditions hold, and we set 

< = inf{j : ^ G [y*, ^)}, i = 1, . . . M 7 . 

Denote by r| the /c-th return of /3 e to the origin. We also set of := and crf^ +1 := 

~e _ e 
r JV e - T ■ 

Note that each of the variables of +1 — of is bounded below in distribution by rj? ^ 2e p 
and recall that for every c > 0, ^/g] converges in distribution to r c . Since r c > a.s., we 
may choose 5 G (0, k/2) so small that, for e small, 

P[<r e i+1 - at > 26 for every i G {1, . . . , M 7 }] > 1 - |. (4.11). 

Write for the intersection of the set i£ e with the event considered in (4.11). Notice 
that on E' £ we can choose a finite sequence = Sq < sf < ■ • • < s £ Ke = r e in such a way 
that 5 < Sj — s|_! < 25 < «, for every j G {1, . . . , _ftT e }, and each interval [s^^Sj) is 
contained in exactly one interval [crf^, a k)- 

We use the sequence (sf ) to get an upper bound on 6{e, 6) on the event E' e . First 
observe that for j G {1, . . . , K e }, 

sup d{W^W £ s> )< sup sup sup|W s e (t)- W s e ,(t)|. 

s , s '6[^_ 1)S p s,s'e[^_ 1)S p s.s'el^.i.sp *>o 

The first term on the right hand side is bounded above by p < rj/h by the definition of E e 
(cf (4.10)) and the property s e - — s e j_ 1 < k. To bound the second term, let s, s' G [sj_i, s|) 
and consider first the case when 

m e (s,s') := inf ^ > 0. 
Then W s e (t) = W s e ,(t) for every t G [0, m e (s, s')], and thus 

sup \wi{t) -wi,{t)\ 



t>o 



< sup \W e s {t) - W*(m e (s,s'))\ + sup \W e s ,{t) - W e s ,{m £ {s, s'))\ 

m s {s,s')<t<(3 e s rn s {s,s')<t<f3 £ a , 

< 2 A 

~ 5 
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again by the definition of E e (cf (4.9) and (4.10)). The case m £ (s,s') = is analogous, 
but we now get the additional term |W/(0) — W £ , (0)|. However, by construction, s and 
s' belong to the same interval [cr|_ 1; a|) and thus W £ (0) and W £ ,(0) belong to the same 
[yk, Zk), which implies that | W/(0) - W £ s , (0) | < 7 = i]/5. Finally, for every j e {1, . . . , K e }, 
we get the bound 

sup d(WZ,WZ,)<^<ri 

s,a'e[s|_i,s|) 

on It follows that, for £ small, 

P[9(e,5)>r ] ]<P[(E' £ y]<^<r 1 . 

This completes the proof of (4.8) and of the tightness of the sequence W e . 

The remaining assertions of Proposition 4.2 are easy, (i) is clear since f3 must be the 
weak limit of f3 e . (ii) follows from the analogous property for W e , and a similar argument 
applies to (iii). □ 

4.3 Tightness of the reflected historical processes 

Recall that the historical process for the ^-reflected system is the process with values 
in M f (W) defined by 




It is easy to verify that Y £ has right-continuous paths with left limits. The following 
theorem is a slightly more precise version of Theorem 1.1. 

Theorem 4.3. The sequence of the laws C Y ofY s , e e S, is tight in the space of probabil- 
ity measures on D([0, 00), M/(W)) and any limit law is supported on C([0, 00), M/(W)). 
Suppose that Ly is the limit of a subsequence of C Y - By passing to a further subsequence 
of e 's, if necessary, we may assume that the laws £fy of W s converge to a law Cw ■ Then 
one can construct on some probability space processes Y and W with distributions Ly and 
Cw , resp., related by 

Y t = £ dl\ 5~ s , 

where (L* , t > 0, s > 0) denote the local times of the process j3 s := ( w , and r = inf{s > 
0:L° = a}. 

Proof. By Proposition 4.2, the laws of W s , e G £ are tight. Hence, from any subsequence 
of £, we can extract a further subsequence Sq along which W e converges in distribution. 
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We can in fact obtain more. For every e > and t > 0, denote by Tf, Tf the random 
measures on R + defined by 

(rf,y>> = r dL^^s) , (n,^ = r dL*/y{ S ). 

Jo Jo 
(We have r e = r e but we prefer to keep a different notation here.) Also define T t by: 

(T u <p) = [ dL*Ms). 
Jo 

As a consequence of Lemma 2.1, we know that 

e— >0 

uniformly on [0, oo) 2 , a.s. If we replace the pair (/3 £ Ar e, T e ) by (f3 s Af£ , T e ) this convergence 
still holds in distribution in C(R+,R) x D(R + , M/(R+)). From this observation and 
standard arguments, we have the joint convergence 

(W s ,p s ,T s ) (W,p,T) (4.12) 

where 

(f u <p)= [ d&Ms), 

Jo 

with the notation introduced in the theorem. 

By the Skorohod representation theorem, we can replace for every e G £q the triplet 
(W s , (3 e ,T e ) by a new triplet having the same distribution, in such a way that the con- 
vergence (4.12) now holds almost surely. Without risk of confusion, we keep the same 
notation for the new triplets. We claim that we have then 

Yf = I fUds) 5~ E — — > f f t (ds) 5~ = Y t (4.13) 

uniformly on compact subsets of R+, a.s. Clearly Theorem 4.3 follows from (4.13) and the 
fact that the limiting process Y that appears in (4.13) is continuous. Both (4.13) and the 
latter fact are immediate consequences of the convergence (4.12) (now assumed to hold 
a.s.) and the following "elementary" lemma, whose proof is left to the reader. 

Lemma 4.4. Let (7", n E N) be a sequence in D(R + , M/(R + )). Assume that 7™ converges 
as n — > 00 to 74, uniformly on every compact of~R + , that t — > 7 t is continuous and that 
the measure "ft is diffuse, for every t G R+. Let E be a Polish space and let (f n , n G N) 
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be a sequence in D(R + , E) that converges to f in D(R + , E). For every integer n G N and 
every t G K+, Ze£ z/ t n G Mf(E) be defined by 

Then converges as n — > oo ; uniformly on compact subsets of~R + , to the measure u t 
defined by 

vt = J lt{ds) 5 f(s) . 

Furthermore, the mapping t — > v t is continuous. □ 

Remark. We do not know whether the limit law of the sequence C Y in Theorem 4.3 is 
unique. A positive answer would give the convergence in distribution of the processes Y e . 
We can also formulate the problem in terms of the reflected snake. Is there a unique (in 
law) process W satisfying properties (i) - (iii) of Proposition 4.2 and such that 

dZt ° S w s (t) 

is a super-Brownian motion started at fx ? 

5. Path properties of the reflected historical process 
5.1 Preliminaries 

Throughout this section, we consider a process Y which is a weak limit of the processes 
Y e as e — > 0. According to Theorem 4.3, we may and will assume that Y is constructed 
together with the reflected Brownian snake W, in such a way that, for every t > 0, 

Y t = f Q dL\8~ s 

where (L*, t > 0, s > 0) denote the local times of the process l3 s := , which is (twice) a 
reflected Brownian motion stopped at time r — mf{s > : L° s = a}. 
The process 

is the weak limit of the processes X s = X s and therefore must be a super-Brownian motion 
started at [i. 

Let us recall the two key properties of the reflected snake W (cf Proposition 4.2): 
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• Monotonicity property: Almost surely for every s < s' we have W s (t) < W s >(t) for 
every t G [0, fi s A(3 S /]. 

• Snake property: The set of discontinuities of the mapping s — > W s is contained in 
the zero set of (3. Furthermore, if s < s' belong to the same connected component of the 
complement of the zero set, we have 

W s (t) = W s ,{t) for every t G [0, inf r ] . 

r€[s,s'] 

In order to state a useful preliminary result, we introduce some notation. Let us fix 
t > 0, and denote by (aj, bf), i G It the excursion intervals of (3 above level t (equivalently, 
these are the connected components of the open set {s > : (3 S > t}). Note that the index 
set It may be empty. For each % G It, denote by e\ the corresponding excursion 

e*( S ) = P(a*+8)Ab* -*» s > °- 

By the snake property of W, we have 

W s {t) = WAt)=:zl VaG[a*,6S]. 

We denote by n(de) the Ito measure of positive Brownian excursions. We normalize 
the measure n(de) by declaring that the Poisson point process of excursions from 0, i.e., 
the family of points (L° , e°), has intensity dsn(de). 

i 

Proposition 5.1. Conditionally on X t , the point measure 

is Poisson with intensity X t (dz)n(de). Consequently, for every Borel subset A o/R, the 
process 

r^Zy = J Y t+r (dw) l A (w(t)) 

is a Feller diffusion started at X t (A). 

We recall that the Feller diffusion is a diffusion process Z on R + whose transition 
kernels are characterized by the Laplace transform: i?[exp(— XZ t )\Z = z] = ex.-p(—zu t (X)) 
where 

»<< A > = l4~pt' 

The total mass process {X t , 1) = L% is a Feller diffusion started at a. 
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Proof. We denote by 77 the right-continuous inverse of the function r — > L* . Note that 
Tr < oo iff r < L~ = (X t , 1). We can rewrite the definition of X t as 

(X t ,¥>> = / ^r^oft)). (5.1) 
Jo 

We also set for every r > 0, 

and we let 7r be the right-continuous inverse of the function r — > . Finally we set 

=(3 w - t, for every r G [0, 

We then claim that, conditionally on {L~ = x}, the process , < r < is 
a reflected Brownian motion started at and killed at the first hitting time of x by its 
local time at level 0, and is independent of the process (W. w (t),0 < r < L%). Except for 
the independence statement, this is a familiar property of linear Brownian motion: See 
e.g. Section VI. 2 of [RY]. To get the independence property, observe that the analogue of 
the process (3^ for the ^-reflected system codes (in the sense of Section 2) the genealogy 
of the descendants of particles at time t. On the other hand, if t 6 ^ 1 ' denotes the right- 
continuous inverse of L £,t , the process (W s £ (t) (t), r > 0) just enumerates in increasing 
order the positions of the particles alive at t. The required independence is thus clear at 
the discrete level of the £-reflected system, and it is preserved under the passage to the 
limit (4.12). 

To complete the proof, write if" 1 for the local time at of (3^> at the beginning, or 
the end, of excursion e\. Note that r ( t) = b\ and thus 

z®=W rW (f). (5.2) 

/<*> 

The point measure ^<W*) n * s the excursion process of the process /3^\ Hence, condi- 

tionally on {L\ = x} : this point measure is Poisson with intensity l[ 0jX )(£)d£n(de) and is 
independent of (W.(t) (t), < r < L|). The first part of the lemma then follows from this 
property, (5.2) and (5.1) (which just says that X t is the image of the measure l r j( t )Al)dl 

under the mapping i^W (*)(£)). 

T t 

To get the second assertion of the lemma, note that by the definition of Y t+r , 
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where i r (e\ ) denotes the total local time of excursion e\ at level r. By the first part of 
the proposition, conditionally on X t , the random measure 

is Poisson with intensity X t (A) n(de). Hence, conditionally on {X t (A) = x} : the process 
(Zp A ,r > 0) has the same law as {L r Tx ,r > 0), and the desired result follows from the 
celebrated Ray-Knight theorem on Brownian local time. □ 

Remark. We could easily sharpen the statement of Proposition 5.1 by conditioning on 
Y t , or even on (Y u , u <t) rather than on X t . We will not need these refinements. 

5.2 A priori estimates 

By [KS] or [R], we know that, almost surely for every t > 0, the measure X t has a 
continuous density xt(y) with respect to Lebesgue measure on M, and the family (xt(y),t > 
0,|/ £ 1) is jointly continuous. Some of our results will be proved under the following 
additional assumption: 

Assumption (H). The measure \i has a continuous density xo(y) with respect to Lebesgue 
measure. 

Under (H), the family (x t (y),t > 0, y G K) is jointly continuous (see Theorem 8.3.2 in 
[Da]). 

In order to simplify the statements of the results in this subsection we introduce a 
constant a. All the results hold for a = 0, assuming (H). Without this assumption, the 
results hold for any fixed strictly positive a. 

For every t > 0, r > and z G K, we set 

ipt,t+r{z) = sup{W s (t + r) : (3 S > t + r and W s (t) < z}, 
with the usual convention sup0 = — oo. We also consider the symmetric quantity: 
$t,t+r(z) = ™f{W s (t + r) : [3 S > t + r and W s (t) > z}, 

Proposition 5.2. Let n G (0, -|) and c > 0. Then, almost surely, one can choose Sq > 
small enough so that, for every 5 G (0, do), t > a and z G M, the condition x t (z) > c 
implies 

il>t,t+6(z)>z-6t-i. 
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Proof. For every t > and z G R set 

7*' 2 = inf{s > : & > t and W^t) > z}, 

with the convention inf = r. Using the formula for X t in terms of W, and then the 
monotonicity property, we get 

pT pT 

X t ((-oo,z]) = dLl l { ~ s(t)<2} = / dL\ l {s<7 t, z} = L* t ,». 

t/ */ 

On the other hand, if s < and /3 S > t+8, we have W s (t) < z and W / s (t + 5) < ipt,t+s(z)- 
Therefore, 



X t+ ,((-oo, Vt,t+«(*)]) = jf dL*+* 1 { 



W.(t+«)<^t,t+«(*)} - T*' z - 



Thanks to the Holder continuity of Brownian local time in the time variable, we can choose 
Si > so small that, for every 8 G (0, Si], t > and zGR, 

L$H >L^ Z -S^. 
By combining all these facts we obtain for every S G (0, t > and zGK, 

JW(-<x>,V M+ *(z)]) >*t((-oo, *])-<**-". (5.3) 

Note that the set {(£, y) : x t (y) > 0} is contained in the graph of X and is thus relatively 
compact. By uniform continuity, we can choose S2 > small enough so that, for every 
t > a and z£l, the condition x t (z) > c implies that x t +s(y) > § for all S G [0,82] and 
y G [z — 82, z + 62]- In particular, if < r < 82 and 8 G [0, 82], 

X t+5 ((-oo, z-r])< X t+s ((-oo, z])-^r. (5.4) 

The proof of the following simple estimate for super-Brownian motion is postponed to the 
appendix. 

Lemma 5.3. Almost surely there exists 83 > such that, for every t > a, z G R and 
8e(0,8 3 ), 

\X t+s ((-oo, z]) - X t ((-oo, z])\ < 6$-". (5.5) 



To complete the proof of Proposition 5.2, choose 80 G (0, Si A 82 A 83) and also such 
ioiSe (0,(5 ), 



that |5 2 77 < <^2- Then, if t > a and zGl are such that x t (z) > c, (5.3) and (5.5) give 



X t+5 ((-oo,ip t , t+5 (z)]) > X t+5 ((-oo,z]) -IS" \ 
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Using (5.4) with r = ^6* v , we get 

4 i_ 

X t+5 ((-oo,ip t , t+5 (z)}) > X t+s ((-oo,z- -6* 

which implies 

By replacing r\ with r( G (0, rf) we can get rid of the factor |. □ 

We can immediately use Proposition 5.2 to derive some useful results on continuity 
properties of the paths W s . Note that, if s G (0, r) is such that j3 s > t + r and W s (t) > z, 
the monotonicity property of the reflected snake implies that W s (t + r) > ifj tjt+r (z). Using 
Proposition 5.2 and the symmetric result for if; tjt+r (z), we get the following corollary. 
Recall that we take a = if (H) is assumed to hold and a > otherwise. 

Corollary 5.4. Let r\ £ (0, |) and c > 0. Then almost surely we can choose Sq small 
enough so that, for every t > a and every s G (0, r) such that (3 S > t and x t (W s {t)) > c, 
we have for every r G [t, (t + So) A [3 S ], 

\W s (r)-W s (t)\ < {r-t)^~\ 

5.3 The key technical lemma 

Our aim is to refine the a priori estimates that were derived in the previous subsection. 
To this end, we will need a crucial technical lemma (Lemma 5.7 below), whose proof 
requires coming back to the approximating branching particle systems. Recall the notation 
(W e : f3 £ , Y e ) of the previous sections. A much simplified version of the arguments of Section 
4 yields the convergence in distribution 

(W S ,Y S )^-+(W,Y), 

where W is a minor modification of the Brownian snake of [L2] (to be precise, W is 
obtained by concatenating a Poisson point process of Brownian snake excursions with 
intensity J fi(dy) N y , in the notation of [L2]) and Y is the historical super-Brownian motion 
connected to W via the formula 

Y t = [ dLl8 Ws , 
Jo 

where (L*,t > 0, s > 0) are the local times of the lifetime process f3 s = (w s , which is a 
reflected Brownian motion stopped at time r = inf{s > : L° s = a}. (Our notation is 
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slightly inconsistent with the previous sections, where (3 was not stopped, but this should 
cause no confusion.) 

On the other hand (cf the proof of Theorem 4.3), we may and will assume that there 
is a sequence Sq of values of e such that 

(w £ ,y £ ) J£U (iv,Y). 

By a compactness argument, and replacing the sequence So by a subsequence if necessary, 
we have also 

(W £ ,Y £ ,W £ ,Y £ ) — (W,Y,W,Y). 

By the Skorohod representation theorem, we can for every e G So find a 4-tuple which has 
the same distribution as (W £ , Y £ , W £ , Y £ ) (and for which we keep the same notation), in 
such a way that the previous convergence now holds a.s.: 

(W £ ,Y £ ,W £ ,Y £ ) (W,Y,W,Y). (5.6) 

From now on we will restrict our attention to values of e in the sequence So and assume 
that (5.6) holds. From the equality X £ = X £ , we also have 



J dL f s 5 Ws{t ) = J dL f s S~ s{t) = X t 



and we see that r coincides with r. 

We introduce the following more restrictive version of Assumption (H): 

Assumption (H'). The measure fx has a continuous density xo(y), which is Holder con- 
tinuous with exponent \ — 5, for every 5 > 0. 

As in the case of Assumption (H), in order to be able to use a single statement for a result 
with or without Assumption (H'), we take a = if (H') holds and otherwise we let a be a 
fixed strictly positive constant. We also fix a constant c G (0, 1). 

Let 77,77' ,p be three positive constants, with < n < n' < 1/4 and p G (0, |). For 
every 8 G (0, 1), we denote by E(5) the event on which the following three conditions hold. 

A. For every s > 0, t G [0, a ], and r G [t, (t + 5) A /3 S ], 

\ Ws (r)-W s (t)\<±(r-t)^. 

B. For every t > a and s > such that l3 s > t and x t (W s (t)) > c, we have for every 
r G [t,(t + 6)A(3 a ], _ _ 

\W s (r)-W s (t)\ < {r-t) 1 ^. 
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C. For every t > a, z E R, and y E [z — 5^ v ' 7 z + 5? r} '\, 

\x t (z) - x t (y)\ <\z- y\^~ p . 

Note that the sets E(S) are decreasing in S. We have P[\J n E(2~ n )] = 1. The fact that 
properties A and B hold for 6 small enough follows from the Holder continuity properties 
of the Brownian snake paths (cf (2.7)) and Corollary 5.4 respectively. For property C, see 
Theorem 8.3.2 in [Da] when a > 0. When a = (then (FT) is in force), the desired Holder 
continuity of the densities is easily obtained from formula (8.3.5b) of [Da] by using the 
techniques of [KS]. 

Throughout this subsection, we fix 8 E (0, 1), t > a and We plan to improve 

the estimates obtained on ipt,t+5 (z) m the previous subsection. We set 

7 = <yW 

and we assume that 5 has been chosen small enough so that 7 > 4 5^~ v . Then, for every 
r E [t, t + 8], we set 

X* = J Y r (dw) l{ w ( t )e(z- 1 ,z+ 1 )} $w(r)- 

The random measure X* corresponds, for the historical super-Brownian motion Y, to the 
contribution of those particles alive at time r whose ancestor at time t lies in the interval 
(z — 7, z + 7). Note that X t * is simply the restriction of X t to (z — 7, z + 7). 

Our goal is to compare X^, 5 ((— 00, i^t,t+s(z)]) to 00, z]) in the same way as 

we compared X t +s((— 00, ipt,t+s(z)]) to X t ((—oo,z]) in (5.3) above. Unfortunately, the 
argument has to be significantly more complicated. 

We set for every e > 0, 

r t ,t + s(z) = snp{W s s (t + S) : g > t + 5 and W s s (t) < z}, 

which represents for the £-reflected system the right-most position among those particles 
alive at time t + 5 which are descendants of the particles located to the left of z at time t. 

Lemma 5.5. We have 

ipt,t+d(z) = lim^lt+d( z ) a - s - 

Proof. This is basically a consequence of the convergence of W £ towards W, which entails 
the convergence of f3 £ to f3. We also use the fact that in the definition of ip t ,t+s(z) : i.e., 

ipt,t+s(z) = sup{W s (t + 5) : P s > t + 5 and W s (t) < z}, 
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we can replace the weak inequality j3 s > t + 5 by a strict one, and/or the strict inequality 
W s (t) < z by a weak one. To justify this, note that: 

(a) Almost surely, every s such that f3 s = t + 6 is the limit of a sequence s n such that 
(3 Sn > t + 5 (simply because t + 8 cannot be a local maximum of (3). 

(b) With probability 1, there is no value of s such that W a (t) = z and (3 S > t + 8 (this 
immediately follows from Proposition 5.1). 

We leave details to the reader. □ 

We now introduce a different approximation of ipt,t+s(z)- We consider in the (non- 
reflected) e-system those particles which are located at time t in the interval (2 — 7,2-1-7), 
and the descendants of these particles after time t. With this branching particle system 
(evolving over the time interval [t, 00)), we can associate a reflected system in the way 
explained in Subsection 3.1. We denote by ip^ +g (z) the position in this new reflected 
system of the right-most particle at time t+8, among those particles which are descendants 
of the particles located to the left of z at time t. 

For every r > 0, we set x(t, z, r) = inf {xt(y) : \y — z\ < r} and x(t, z, r) = sup{xt(y) : 
\y — z\ < r}. 

Lemma 5.6. We have 

P[(hmsup{^' t %(2) ^ il>Z,t+6(*)})nE(S) n {x(t,2,^/ 2 ) > c}j < 2exp(-2c5- 1 / 2 ). 



Proof. We introduce the following events: 

A + = {3s > : p s > t + 5 and W s (t) G (z - 5 1/2 , z)}, 

and 

A" = {3s > : P s > t + 6 and W s (t) e(z,z + S 1/2 )}. 
We first verify that a.s., 

((limsup{^ e +5 (2) ± il>t tt+s {z)}) n E(8) n {*(*,*, J 1 ' 2 ) > c}) C (A+ n A") c . (5.7) 

Suppose that A+ n A~ n E(8) n {x(t, z, 5 1/2 ) > c} holds. Then, there exists s 1 > such 
that P S1 > t + 5 and W Sl (t) E (z — 5 1 / 2 , z). From property B in the definition of E(5) we 
also have |W Sl (r) — z\ < 25^~ v for every r G [t,t + 5]. Similarly, there exists S2 > such 
that (3 S2 > t + 5, W S2 (t) G (z,z + 5 1 / 2 ) and \W S2 (r) - z\ < 25^ for every r G [t, t + 5]. 
By the convergence (5.6), the same properties hold for e > small enough, if we replace 
W Si and (3 Si by M 7 ^. and (3 £ s . respectively. 



34 



On the other hand, by property A of the definition of E(5) and the convergence (5.6), 
we have also for e small enough, for every s such that f3 e > t and every r G [£, (t + 5) A 

|W?(r)-W7(f)| 

In particular, if s is such that f3 £ s > t and |W s £ (t) — z\ > 7 > 4 5^ ~ v , we have for every 
re[t, (t + 6)Afi], 

\W e s {r)-z\ >2 5^~ 11 . 

We have shown that, on the event A + fl A - D E(8) fl {x(t, z, 5 1 / 2 ) > c}, provided that 
e is small enough: 

• There exist s\ and S2 such that j3 Sl > t + 5, /3 S2 > t + 5 and 

(t) e(z- z) , \W s Sl (r) - z| < 25*-", Vr e[t,t + 5] 
W s S2 (t) E (z.z + 5 1 / 2 ) , \W s S2 (r)-z\ < 253"", Vr e + 

• For every s > such that f3 e >t and |WJ (t) — z| > 7, 

\W e s (r)-z\ >2 8*~ r> , VrG [*,(* + <J) A . 

These properties allow us to apply Lemma 3.1. In the context of that lemma, the original 
system is the e-system considered after time t, the new (restricted) system consists of the 
descendants of the particles which are located at time t in the interval (z — 7, z + 7), and 
we take I = (z — 25^~ r >, z + 25^~ r >). Lemma 3.1 and the previous properties imply that 
the restrictions of the paths Wg and W £ 2 to [t, t + S] still appear as restrictions of reflected 
historical paths in the new system. Note that in the definition of V^t+a ( z ) > respectively of 
'4 ) t't+s^ z )i we ma y res trict our attention to those reflected historical paths between times t 
and t + S in the original system, resp. in the new system, whose value at time t lies in the 
interval [W/ (t) , z) (this is so because of the monotonicity property of reflected historical 
paths). Any such path is bounded below and above by W% and W% respectively, on the 
time interval [t, t + 5]. By Lemma 3.1 again, the class of paths that we need to consider is 
exactly the same for both the original system and the new one. This is enough to conclude 
that il)t't +5 (z) = ^ £ t+s (z), and we get our claim (5.7). 

It follows from (5.7) that the probability considered in the lemma is bounded above 

by 

P[(A+ n A") c n {x(t, z, s 1/2 ) > c}}. 

By the construction of Y, we have 
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Hence the event A + certainly holds if 

J Y t+S (dw) l{ w (t)e(z-8^/2,z)} > 0- 

It follows that 

P[(A+) c n{x(t,^5 1 / 2 ) > c}} < p[{ J Y t+S (dw) i {t0(t)e(a _ tf i/ a> , )} = oynixfoz^ 1 ' 2 ) > c}], 

and a similar bound holds if we replace A+ by A - . By Proposition 5.1, the last quantity is 
bounded above by the probability that a Feller diffusion started at cb 1 ! 2 vanishes at time 
5. This probability is equal to exp(— 2c5 -1 / 2 ), which completes the proof. □ 

We can now state the key lemma. We fix still another constant rj" G (77', 1/4). 

Lemma 5.7. There exist two positive constants C and k, that depend only on c, 77, rj', rj" 
and p, such that 

p[{\x: +s ((-oo^ t , t+5 (z)]) - x:((-oo,z])\> s^''} 

n E(6) n {x(t, z, 5 1/2 ) > c} n {x(t, z, 7) < c -1 }] < C exp(-<r K ). 

Proof. For every r G [t, t + S], set 

X*' s = J Y^(dw) l{ w (t) e ( z - 7 ,z +7 )}5 w ( r ), 

which represents the contribution at time r of the descendants (in the non-reflected system) 
of particles which are located in (z — 7, z + 7) at time t. From the convergence of Y £ to Y, 
and the fact that f Y r (dw) l{ w (t)=z±-y} = 0, one can easily show that for every r G [t, t + S], 
the measures X*' 6 converge weakly to X*. In particular, a.s. for every y G R, 

hmX^ s ((-oo,y]) =X; +5 ((-oo,y]). 

From Lemma 5.5 and Lemma 5.6, we get that on the set E(5) fl {x(t, z, b 1 ! 2 > c}, we have 
the convergence 

lmiX ( ^((-oo,^; t %(2)]) = X* +s ((-oo,ipt !t+6 {z)]), 

except possibly on a set of measure at most 2 exp(— 2c5~ 1 / 2 ). 

However, by the definition of ty* t ^ +s (z), and the monotonicity property of reflected sys- 
tems, the quantity X^ s ((—oo, ty* t ^ +s (z)}) is equal to e times the number of descendants at 
time t+S of the particles present at time t in (2—7, z), for the ^-reflected system constructed 
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over the time interval [t, oo) from the particles present at time t in (z — 7, 2 + 7). Since the 
law of the branching evolution is the same for the reflected system as for the original one, 
we see that conditionally on {X^' e ((— 00, z]) = ek}, the variable X*^ 5 ((— 00, ip^f +s (z)]) is 
distributed as eZg 6,h \ where Z^ £ ' k ^ denotes a Galton- Watson process with critical binary 
branching at rate e~ x and initial value k. Recall that X^' £ ((— 00, z]) converges a.s. to 
X*((— 00, z]). By standard limit theorems for Galton- Watson processes, 

(x;' £ ((-oo, 2 ]),x;4f 5 ((-oo,^; t %(^)])) (x;((-oc,z]),u), 

where conditionally on X£((— 00, z]) = u, the variable U is distributed as the value at time 
5 of a Feller diffusion started at u. 

Note that X t *((— 00, z]) = -X"t((^ — 7,^]) and that on the set {x(t, 2,7) < c -1 } we 
have ((— 00, z]) < c _1 7 = c~ 1 8^~ ri ' . Elementary estimates on the Feller diffusion, using 
only the form of the Laplace transform of the semigroup (see the appendix for very similar 
estimates) show that 

P[{X;((-oo,z]) < c- 1 **-"'} n -X t *((-oo,z])| > 8^"}} < C exp(-5- K '), 

where the constants C and k' > depend only on c, 77' and rj". 
To complete the proof of the lemma, we write 

P[{\X; +5 ((-oo^ t , t+s (z)]) -Xt((-oo,z})\ > 

n E(5) n {x(t, z, 5 1/2 ) > c] n 2, 7) < c -1 }] 

< 2 exp(-2ctf-V 2 ) + P[{lW^ > tfi"""} 

< 2 exp(-2a5- 1 /2) + lirrnnfP[{|X;4 e 5 ((-oo,^%W]) - X t *((-oo,z])| > 

n{i;((-oo )2 ])<c-^-"'}] 

< 2 exp(-2c5- 1 / 2 ) + P[{X*((-oo, z]) < a" 1 ^'} n {|£7 - X t *((-oo, z])| > <^ _J? "}] 

< 2 exp(-2c5- 1 / 2 ) + C exp(-<r K '). 

□ 

5.4 The main result 

We keep the notation introduced in the previous subsection. The reals t > a, 6 G (0, 1) 
and z G R are fixed for the moment. 

Lemma 5.8. Assume that rj" > |?/ + \p. There exist two constants C and k > 0, that 
depend only on c, rj', r}" and p, such that 

P[{X* +s ((-oc, z\) > X*((-oc, z]) + Si-*"} n E(8) n {x(t, z, 7) < c" 1 }] < Cexp(-<T*). 
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The proof of this lemma is an application of standard techniques in the theory of 
super-Brownian motion. See the appendix for a detailed argument. 

Proposition 5.9. Under the assumptions of Lemma 5.8, there exist two constants Co and 
Ko > 0, that depend only on c, rj, rj', rj" and p, such that 

P[{iP t , t+s (z)<z-- c 5^"}nE(5) 

n {c < x(t, z, 7) < x(t, z, 7) < c" 1 } n {x(t + 5, z, 7) > c}] < C exp(-5~ K °). 

Proof. Our argument is very similar to the proof of Proposition 5.2. We will assume that 
the event E(5) fl {c < x(t, z, 7) < x(t, z, 7) < c -1 } holds. By Lemmas 5.7 and 5.8, we have 
on this set 

X? +s ((-oo, z]) < X*((-oo, z\) + 8-^" < XZ +s ((-oo, lH, t+s (z)]) + 2 5^" (5.8) 

except possibly on a set of probability at most C exp(— S~ K ) + C exp(— 5~ K ). 

On the other hand, condition A in the definition of E(5) (and the fact that 7 > 4 S^~ v ) 
ensures that the measures X^ +5 and X t+ $ coincide over the interval (z — z+ Hence, 
on the event {x(t + 5, z, 7) > c}, we get 

X; +s ((-oo, z})-2 8^" > X: +s ((-oo, z - ^ W]), 

provided that 5 is small enough so that ^Si~ v < ^. On the set where (5.8) holds, we get 

x; +5 ((-oo^ M+ ^)])>x; +5 ((-oo,z-^t-^]), 

and the desired result follows. □ 

We now come to the main result of this section, which is a refinement of Corollary 
5.4. Recall our conventions concerning a — this constant is equal if (H') is assumed to 
hold and otherwise a is a fixed strictly positive constant. 

Theorem 5.10. Let A > and c G (0, 1). Then a.s. we can choose Sq small enough so 
that, for every t > a and every s G (0, r) such that j3 s > t and x t (W s (t)) > c, we have for 
every r G [t, (t + S ) A (3 S ], 

\W a (r)-W a (t)\ < (r-t)i~ x . 
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Proof. We can choose 77, ?/, rj" with < 1] < rj' < rj" < A and p G (0, \) such that 
the assumptions of Lemma 5.8 hold. We then apply the estimate of Proposition 5.9 with 
5 = 2~ n (n large enough) to all reals t G [a,n], z G [— n, n] of the form t = k2~ n , 
z = j2~ n . We have already observed that P[\J n E(2~ n )] = 1. Furthermore, if we assume 
that c < x t (z) < c" 1 we will have x(t, z, 2" n ^-V)) > c / 2? x(t, z, 2" n (W)) < 2 / c , and 
x(t + 2 _n , z, 2~ n( -^~ v )) > c/2, for all n sufficiently large (depending on u; but not on t and 
z). Then, by combining the estimate of Proposition 5.9 with the Borel-Cantelli lemma, we 
obtain the following property: There exists an integer no(u) such that for every n > no(u), 
for every t = /c2 _n , z = j2~ n with t G [a,n], z G [— n, n], the condition c < xt(z) < c _1 
implies 

^t,t + 2-n{z)>z-{2- n )i- x . 

Since the densities x r (y) are bounded over [a, 00) x R, a simple argument shows that we 
can drop the condition xt(z) < c _1 in the previous assertion. 

Then, if s > is such that (3 a > t + 2 _n , where t is of the form t = /c2 _n , we let 
z = j2~ n be such that z < W 8 (t) < z + 2~ n . If n is large enough (again independently of 
the choice of s and t), the condition Xt(W s (t)) > 2c will imply Xt(z) > c. Then, by the 
definition of tfj t ,t+s(z) and the preceding estimate, 

W s (t + 2~ n ) > («) > W s (t) - 2~ n - {2~ n )i-\ 

Thanks to this observation and a symmetry argument, we obtain that a.s. for n large 
enough, for every t > a of the form t = k2~ n and every s > such that j3 s > t + 2~ n and 
x t (W s (t)) > 2c, 

\W 8 (t + 2~ n ) -W s (t)\<2 (2" n )|- A . 

The statement of Theorem 5.10 now follows easily thanks to the usual chaining argument. 

□ 

Theorem 1.2 is an immediate consequence of Theorem 5.10. Note that, by the repre- 
sentation formula for Y in terms of W, the set supply is contained in {W s ;f3 s = t}, for 
every t > 0, a.s. The comments following the statement of Theorem 1.2 are justified by 
Proposition 5.1. 

6. Branching points 

In this last section, we prove Theorem 1.3. As in Section 5, we assume that the process 
Y is constructed together with the reflected Brownian snake W, in such a way that we 
have the representation formula 
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We need a preliminary lemma. If si < S2, we set m(si, s 2 ) = inf s6 [ Sl S2 ] j3 s . 

Lemma 6.1. Almost surely, for any t > and any s\ < S2 such that j3 Sl = (3 S2 = t and 
< m(si, S2) < t, we have 

(W ai (m(si,s 2 ))) >0. 



Proof. Let a > and let A > 1 be an integer. Write Ea for the event Ea = {Q C 
[0, A] x [—A, A]}, where Q is as above the graph of X. It is enough to prove that a.s. on 
Ea, the following holds: 

(P) For any t > a and s± < S2 such that j3 Sl = (3 S2 = t and < m(si, S2) < t — a, we have 
(W ai (m(si,s 2 ))) >0. 

We first introduce some notation. Let e be an excursion, that is a continuous function 
e : R+ — > R+ such that e(s) > iff < s < a(e), for some a(e) > 0. Set 



and, if T a (e) < 00, 



T a (e) = inf{s > : e(s) = a} 



L a (e) = sup{t > : e(t) = a}, 
MJe) = inf e(s) 

T a (e)<s<L a (e) 



By convention we take M a (e) = if T a {e) = 00. 

Let r > 0. Recall the notation I r and e\, z\, i E I r introduced before Proposition 5.1, 
and for every c > and 5 G (0, a), set 

N*(a,c) = ^ 1 {x r {z r i )<c} 1 {0<M a (ep<<5}- 

i&Ir 



Proposition 5.1 allows us to conclude that, 
E[N 5 r {a,c) 1 EA ]<E[J21 



iei r 



{k r |<^} 1 {x r (zO<c} 1 {0<M a (e r i )<S} 



= E\ I dzx t (z) l {XrW < c} n(0 < M a (e) < 5) 



< 2cA5a~ 2 , 

using the easy formula n(0 < M a (e) < S) = 5a~ 2 . We apply this estimate with S = 1/k 
(k large enough) and r = j/k for all j = 1, 2, . . . , Ak. It follows that 



E 



i=i 



< 2c A 2 a -2 . 
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In particular, if E k (a, c, A) denotes the event {3j > 1 : N^j^(a, c) > 1} fl Ea, we have 

phiminf E k (a,c,A)] <2cA 2 a~ 2 . (6.1) 

Suppose that property (P) fails. Then, we may find t > a and si < s 2 such that f3 Sl = 
(3 S2 = t and < m(si,s 2 ) < t — a, and furthermore ^ m ( SljS2 ) (W Sl (m( Sl ,s 2 ))) = 0. We 
take j such that j/k < m(si, s 2 ) < (j + l)/k, and observe that Xj/ k (W Sl (j /k)) < c for all 
k sufficiently large, by the joint continuity of densities. Hence by considering the excursion 
of (3 above level j/k that contains s\, we see that N^j^(a, c) > 1 for all k large. Therefore, 
if A) denotes the event on which (P) fails, we have 

P[F(a, A) n E A ] < P[liminf E k {a,c, A))] < 2c A 2 a -2 . 

Since c was arbitrary, we have P[F(a, A) fl Ea] = 0, which completes the proof. □ 
Proof of Theorem 1.3. The representation formula for Y t implies that 

suppY t = {W s :p s = t}. 

(Note that the set on the right hand side is closed, by the continuity properties of W.) 
Hence if w\ and 102 belong to suppYj and wi 7^ W2, we can find si and s 2 such that 
(3 Sl = P S2 = t, and W Sl = wi, W S2 = W2. With no loss of generality, we can assume 
si < s 2 . We claim that 

m( Sl ,s 2 ) = inf{r G [0,t] : Wl (r) ^ w 2 (r)}. (6.2) 

The inequality m(si,s 2 ) < inf{r G [0, t] : wi(r) 7^ w 2 (r)} is immediate from the snake 
property (when m(si, s 2 ) = there is nothing to prove). On the other hand, if we assume 
that there is a rational r G (m(si, s 2 ), t) such that W Sl (r) = W S2 (r), then the monotonicity 
property implies W s (r) = W Sl (r) for every s G [s\, s 2 ] such that j3 s > r. Hence, 

X r = / dL^~ > / dL^~ = (L r S2 - L r Si )8 Wai{r) , 

JO J si 

which gives a contradiction since L r S2 — L r Si > by standard properties of linear Brownian 
motion. 

From now on, we assume m(si,s 2 ) > 0. Note that we have also m(si,s 2 ) < t since 
we assumed that wi 7^ w 2 . By Lemma 6.1, we have x m ^ SliS2 ){W Sl (m(si, s 2 ))) > 0. By 
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monotonicity (and the fact that the measure X r gives no mass to singletons), we get for 
every r G (m(si, s 2 ),t), 



s 2 _ rW 32 (r) 

dL r s = X r ((W Sl (r),W S2 (r)))= I dzx r (z), 

si JW sl (r) 

and by the continuity of densities, it follows that 

W (r) - W (r) ~ 
, ljm S2 y sA ' =x m(si , S2) (W Sl (m(s 1 ,s 2 )))>0. (6.3) 

■52 * 1 

Thanks to (6.3), the behavior of w\(r) — w 2 (r) as r j 7 Wl , W2 = S2) is reduced to that 

of L^ 2 — L^. Write sq for the (unique) time in (s\, s 2 ) such that (3 So = m(si, s 2 )- Standard 
results on Brownian path decompositions show that, for events that depend only on the 
asymptotic cx-field at time 0, the processes {(3 So - u — (3 So ,u G [0, so — si]} and {(3 S()+U — 
f3 So: u G [0, s 2 — «o]} behave as two independent 3-dimensional Bessel processes. It follows 
from this and the Ray-Knight theorem that the process 5 — > £^( Sl,S2 )+ 5 _ £^( s i- s 2)+<5 
the same local path properties (for 5 close to 0) as the sum of two independent squares of 
2-dimensional Bessel processes, which is the square of a 4-dimensional Bessel process. If 
5 — > R$ is the square of a 4-dimensional Bessel process, the law of the iterated logarithm 
shows that 

r Rs 1 
limsup — = 1. 

sio 2<5 log I logd| 

On the other hand, from the well-known rate of escape for Brownian motion in space 
(Theorem 6 in [DE] combined with time-inversion), we have for a > 0, 

r Rs 

lim — — = 00. 

sio S\ log^l" 1 " 

We have just argued that the same properties hold if we replace Rs with £^( Sl ' S2 ) +<5 _ 
£m( Sl)S2 )+5 Thig and imply Theorem 1 3 n 

Appendix 

Proof of Lemma 5.3. For a fixed value of z, the estimate of Lemma 5.3 follows from 
[P2]. As we need uniformity in z, we will provide a detailed argument. Recall the notation 
from Subsection 5.2, and especially the conventions concerning the constant a. Recall that 
Q denotes the graph of X and for every integer A > 1 consider the event 

E A = {QC [0, A] x [-A, A}; sup x t (y) < A}. 

t>a,yEM 
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Note that P[Ea] | 1 as A f oo. (We use assumption (H) when a = 0.) The key step of 
the proof is to show the following inequality for alH > a and zel, 



P[{\X t+d ((-oo, z}) - X t ((-oo, z})\ > 6^} HE A ]<C exp(-<T K ) 



(^41), 



where the constants C and k > may depend on A but not on t, z and 5. To prove 
(Al), we may apply the Markov property at time t and reduce the problem to the case 
t = 0. More precisely it is enough to consider a super-Brownian motion V = (Tt,t > 0) 
with initial value To(dz) = g(z)dz, with a function g bounded above by A and such that 
f g(z)dz < 2A 2 , and to prove that for every 5 G (0, 1), 



Let us first bound P[r 5 ((-oo, 0]) < r ((-oo, 0]) - 8^]. We know that for every A > 0, 

£[exp(-Ar 5 ((-oo,0]))] =exp(-(r ,u«)), 
where ut(z) solves the integral equation 



if B is a linear Brownian motion started at z under P z . The integral equation gives the 
bound 



P[|r 5 ((-oo, 0]) - r ((-oo, 0])| > 6*-"] < C exp(-cT K ). 



(A2) 



1 r r i 

(z) + -E z / u t - r {B r fdr =XP z [B t <0], 



1 Jo 



ut(z)>XP z [B t <0]-—t. 



We use this bound in the following estimates, 



p[r 5 ((-oo,o])<r ((-oo,o])-^- T '] 

< exp(-A5^- ?? + Ar ((-oo,0])) J E[exp(-Ar 5 ((-oo,0]))] 

<exp(-A^-^ + y5(ro,l))exp(A(r ((-oo,0])- j dzg{z)P z [B 6 < 0])) 



Note that for every e > 0, 




with a constant C £ depending only on e and A. By choosing A = 7 2+ e with < e < 77, we 
arrive at the desired estimate for P[Ts((— 00, 0]) < r ((— 00, 0]) — 5^~ v ]. Slightly different 
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arguments apply to P[Tg((— oo, 0]) > r ((— oo,0]) + 5^ r >]. In fact, it is easier to observe 
that 

P[r 5 ((-oo, o]) > r ((-oo, o]) + s*-"] 

< P[(T 5 , 1) > (To, 1) + + P[T 5 ((0, oo)) < r ((0, oo)) - m 

We have just shown how to bound the second term on the right hand side of (A3). As for 
the first term, we need simply recall that (T t , 1) is a Feller diffusion and use the fact that 
for A G (0, f) 

E[exp(A(r,,l))]=exp(^|il). (A4) 

This immediately leads to the estimate needed to complete the proof of (A2) and (Al). 

From (Al) and the Borel-Cantelli lemma, we get that a.s. there is an integer n (cu) 
such that, for every n > n , for every t > of the form t = j2~ n and every z G K. of the 
form z = k2~ n , we have 

\X t+2 - n ((-oc,z])-X t ((-oc,z})\ < (2~ n )^. 

Note that for every fixed z, the process t — > X t ((— oo, z]) has continuous sample paths a.s. 
(see e.g. Corollary 6 in [P2]). The proof of Lemma 5.3 is easily completed thanks to this 
observation, the preceding bound and the usual chaining argument. □ 

Proof of Lemma 5.8. This is very similar to the proof of (Al) above. Note that the 
process (X t * +r ,0 < r < S) is a super-Brownian motion started at X%, which is simply 
the restriction of X t to [2 — 7,2 + 7]. Thanks to this observation and the definition of 
E(5), we see that it is enough to prove the following statement. Let V = (r r ,r > 0) be 
super-Brownian motion with initial value To(dz) = g(z)dz. Assume that the function g 
vanishes outside [—7,7] and that c < g(z) < c -1 and \g{z) — g{z')\ < \z — z'\^~ p for all 
z, z' G [—7, 7]. Then, 

P[T s ((-oo, 0]) > r ((-oo, 0]) + 6*-""] < Cexpi-S- 15 ), (A5) 

where the constants C and k depend only on c, ?/, rj" and p. 
In a way similar to (A3) we first write 

p[r s ((-oo, 0]) > r ((-oo, 0]) + <y W] 

< P[(T S , 1) > (r , 1) + Ui-^"} + P[r s m 00)) < r ((o, 00)) - itfi-""]. 
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Thanks to (A4), we see that, for A < 2/5, 

P[(T 5 , 1) > (To, 1) + \S-^"] < exp(-A((r , 1) + \s^"))E[e x ^\ 

A s „ ( \ 2 (T ,l}5/2 ^ 
= exp(--5 4 Mexp( x _ j. 



Since (r , 1) < 2c _1 7 = 2c~ 1 55~ ?? ', we get a bound of the desired form by taking A = 
S~i +e with r\" > e > rj'. 

For the other term, we proceed as in the proof of Lemma 5.3: 

P[r 5 ((0,oo)) < r ((0,oo)) - < exp(A(r ((0,oo)) - ±6*-*"))E[e- xr '« ' oo \ 

and £[ e -AiM(o,oo))] = eX p(-(r , u 5 )), with u s (y) > XP y [B 5 > 0] - ±A 2 5. It follows that 

p[r 5 ((o,oo))<r ((o,oo))-^t-V] 



< 



-| \2 . roc p 

eX j>(--\6i-i" + y 5(r , 1)) exp (A( J dzg(z) - J dzg(z)P z [B s > 0])) 



< exp(--\di- v " +c" 1 A 2 5 7 ) exp(4A 7 i- p ), 

where in the last line we used our assumption that \g(z) — g(0)\ < \z\^~ p to bound 
J °° dzg(z) — J dzg(z)P z [Bs > 0]. In view of the assumptions of Lemma 5.8, we can 
now choose A = 5i~ e : with rj" > e > §?/ + |, and we arrive at a bound of the desired 
form. This completes the proof. □ 
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